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Abstract. We study a model of fast magnetic recon- 
nection in the presence of weak turbulence proposed by 
Lazarian and Vishniad (1 1999b using three-dimensional direct 
numerical sim ulations. The mode l has been already success- 
fully tested in iKowal et alJ ( 120091) confirming the dependen- 
cies of the reconnection speed Vrec on the turbulence injec- 
tion power Pinj and the injection scale Imj expressed by a 

constraint Vrec ~ -P./ ,- LJj and no observ ed dependency on 
Ohmic resistivity. In IKowal et al.l 12009), in order to drive 



turbulence, we injected velocity fluctuations in Fourier space 
with frequencies concentrated around kinj — l/linj, as de- 
scribed in Alvelius (1999). In this paper we extend our pre- 
vious studies by comparing fast magnetic reconnection under 
different mechanisms of turbulence injection by introducing 
a new way of turbulence driving. The new method injects 
velocity or magnetic eddies with a specified amplitude and 
scale in random locations directly in real space. We pro- 
vide exact relations between the eddy parameters and tur- 
bulent power and injection scale. We performed simulations 
with new forcing in order to study turbulent power and injec- 
tion scale dependencies. The results show no discrepancy 
between models with two different methods of turbulence 
driving exposing the same scalings in both cases. This is 
in agreement with the Lazarian & Vishniac (1999) predic- 
tions. In addition, we performed a series of models with 
varying viscosity i/. Although Lazarian & Vishniac (1999) 
do not provide any prediction for this dependence, we report 
a weak relation between the reconnection speed with viscos- 
ity, Vrec - V-^^^- 



1 Introduction 

Magnetic fields are observed in many astrophysical 
objects and usually play an important or even crucial 



role in their dynamics (see ICrutcherl Il999[ iBeckl |2002[ 
Valleel 1 19971 Il998l e.g.). They are a key ingredient of 



astrophysical processes such as magneto-rotational insta- 
bility, magnetic dynamo, transport and acceleration of 
cosmic rays, accretion disks, turbulence, solar phenom- 
ena, gar nma r ay bursts, etc. ffialbus and Hawlevl 1998 : 
Parker 
201 IT 



1 19921: iHa nasz et al., 2009; Kulpa-Dvb el et al 
Schlickeiser and Lerchet 1 19851 iMelrosei 120091 
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lElmegreen and Scaloll2004l: iKotera and Ohntoll201 lb 

Magnetic fields are solenoidal and evolve only through 
changes in the curl of the electric field. In the limit of zero re- 
sistivity the topology of the field lines is a constant of motion 
and the magnetic flux threading any fluid element is constant. 
Generating large scale magnetic fields re quires some kin d of 
battery effect, like the Biermann battery (lKhannaill998b and 
generating strong large scale magnetic fields requires a dy- 
namo (see IParken 1 19921 for example). In the limit of very 
small resistivity, which is typical for astrophysical objects, 
the magnetic flux is "frozen in" and magnetic field lines will 
resis t passing throu gh one another or changing their topol- 
ogy ( lMoffatlll978h . Due to the presence of plasma motions, 
in particular turbulence, this would result in a very complex 
tangle of field lines in real objects, with negligible large scale 
magnetic flux. However, observations indicate that the mean 
and turbulent components of magnetic fields in many astro- 
physical objects are of similar strengths (see lBeckll2002l for 
example). This implies the existence of a process which can 
violate the frozen-in condition on dynamical time scales, i.e. 
fast magnetic field reconnection. 

The first analytic model for magnetic reconnection was 
proposed independently by Parker ( 1957|) and|Sweet ( 1958|). 
Sweet-Parker reconnection has the virtue that it reUes on a 
robust and straightforward geometry. Two regions with uni- 
form magnetic fields are separated by thin current sheet. The 
speed of reconnection is given roughly by the resistivity di- 
vided by the sheet thickness. However, the plasma in the cur- 
rent sheet is constrained to move along the local field lines. 
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and is ejected from the edge of the current sheet at the Alfven 
speed, Va- Since the width of the current sheet Umits the flux 
of expeUed plasma, the overall reconnection speed is reduced 
from the Alfven speed by the square root of the Lundquist 
number, S = LVa/tj, where 77 is the resistivity and L is the 
length of the current sheet. In most astrophysical contexts 
S is very large and the Sweet-Parker reconnection speed, 
Vsp ~ VaS^^^^, is negligible. Fast reconnection requires 
that the dependence on rj be erased. Given the simplicity of 
the Sweet-Parker setup, this requires that the simple geome- 
try of the current sheet must be broken. 

The realization that Sweet-Parker reconnection is inade- 
quate to explain magnetic reconnection in an astrophysical 
context was immediately apparent, and gave rise t o decades 
of research on models of fast reconnection (see iBiskamp , 
I2OOO: Priest and Forbesl uOOOl for reviews). The first pro- 
posal was to replace the current sheet with an X-point con- 
figuration, so that the "sheet" thickness and length are com- 
parable. This is the basis for Petschek's model of fast re- 



connection (iPetscheld 119641) . However, a dynamically self- 
consistent X-point requires that the outflow prevent a general 
collapse into a narrow current sheet. Otherwise we would 
expect that the same bulk forces that brought the magnetic 
field lines together would lead to Sweet-Parker reconnec- 
tion. IPetschek (1964) proposed that slow-mode shocks on 
either side of the X-point would serve this purpose. More- 
over, those shocks are responsible for converting most of 
the magnetic energy into the heat and kinetic energy. The 
X-point in this model has an overall size which depends on 
resistivity, but since the magnetic field decrease logarithmi- 
cally when approaching the current sheet (due to the assump- 
tion of the current-free magnetic field in the inflow region), 
the resulting reconnection speed is some fraction of Va- Nu- 
merical simulations with uniform resistivity (Biskamp, 1996) 
have showed that in the MHD limit the shocks fade away 
and the contact region expands into a sheet. The only way 
to make the Petschek configuration stable is by introduc 



ing the l ocal nonuniform resistivity ( Upai arid Tsudal 



Scholer, '1989'; 'Uaai', ' 19921 lYanet al.Ul992[ 



Forbes, 



1977 



l20Q0t 



Shibata and Magara , 20 1^ . 



This leaves the possibility that X-point reconnection is 
stable when the plasma is collisionless. Numerical simula- 
tions dShay et al., 1998, 2004) have been encouraging. How- 
ever, there are several important issues that remain unre- 
solved. First, it is not clear that this kind of fast recon- 
necti on persi sts on scales greater than the ion inertial scale 
(see iBhattac hariee et al., 2003^^_^ Several numerica l stud- 
ies JWang et alt i200k iSmifli et all 120041: iFitzpatrickL |2004 



have found large scale reconnection speeds which depend on 
resistivity, i.e. are not fast. Second, in many circumstances 
the magnetic field geometry does not allow the formation of 
X-point reconnection. For example, a saddle-shaped current 
sheet cannot be spontaneously replaced by an X-point. The 
energy required to do so is comparable to the magnetic en- 
ergy liberated by reconnection, and must be available before- 



hand. Finally, the requirement that reconnection occur in a 
collision less plasma restricts this model to a small fraction 
of astrophysical applications. For example, while reconnec- 
tion in stellar coronae might be describe d in this way, ste llar 
chromospheres are not. More generafly lYamadal (120071) es- 
timated that the scale of the reconnection sheet should not 
exceed about 40 times the electron mean free path. This 
condition is not satisfied in many environments which one 
might naively consider to be collisionless, among them the 
interstellar medium. The conclusion that stellar interiors and 
atmospheres, accretion disks, and the interstellar medium in 
general does not allow fast reconnection is drastic and un- 
palatable. 

An alternative to the X-point geometry is to consider mag- 
netic fields that are chaotic, even if only weakly so. Re- 
quiring the plasma to flow along the local magnetic field 
implies a powerful constraint on re connection only if the 



field lines themselves are laminar. iLazarian and Vishniac 
(1999, hereinafter LV99) proposed a model for fast recon- 
nection which depends on the presence of turbulence, and its 
production of weakly stochastic field lines (also briefly de- 
scribed in Section |2]i. Turbulence is a natural consequence 
of convection in stars and of the magne torotational instabil- 
ity in accretion disks (for a review see iBalbus and Hawlev , 
1998i) . In addition, it is now generally accepted that the 
"Big Power Law in the Sky" indicates the presence of tur- 
bulence on sc ales from te ns of parsecs to thousands of kilo- 
meters (Armstrong etal.l Il995: Chepurnov and LazarianI 
2010h . Among other sources, evidence for this comes 
from studies of atomic hydrogen spec tra in mo l ecular 
clouds and galaxies (see Lazarian andPogosvan , 1999 : 



Stanimirovic and Lazariaru EOOll IPadoan et aL ._ 20061 12009 



Chepurnov et al.l l20ld see also review by LazarianI 12009 
and references therein), as well as r ecent studies of emis- 
sion lines and Farad ay rotation (see iBurkhart et al.i 12010 
Gaensler et al.l I2OI 1 V LV99's model uses the properties of 
turbulence to predict broad outflows from extended current 
sheets. The diffusivity of magnetic field line trajectories in a 
turbulent plasma implies that flows can follow local magnetic 
field lines without being confined to the current sheet. When 
the turbulent diffusivity is less than the ohmic resistivity, this 
model reduces to the Sweet-Parker reconnection model. 
The first test of the LV99 model using three-dimensional 



(3D) MHD simulations was performed in iKowal et al 



(120091) . The main predictions of the model were confirmed. 
In this paper we provide additional numerical evidence of 
magnetic reconnection in turbulent environments by testing 
different mechanisms for injecting turbulence. In §|2] we 
briefly review the LV99 model of reconnection and its the- 
oretical predictions. In ^we describe in details the numeri- 
cal model studied in this paper and new method of turbulence 
driving. Although the initial setup and boundary conditions 
are similar to our pr evious studies, and described in details in 
Kowal et alJ ( 120091) . we briefly describe them here for com- 
pleteness, as well. In ^we present an extensive description 
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of new results obtained from studying our numeric model, 
which we discuss later in ^ In ^we present our main con- 
clusions. 



2 The Lazarian-Vishniac (1999) Model 

The notion that turbulence can influence reconnection rate 
is not unprecedented. The ideas in this regard were dis- 
cussed long before LV99. Howeve r, they f ell sh ort of 
solving the problem. For instance, ISpeiserl (119701) con- 
sidered the effects o f turbu lence on microscopic resistivity, 
Jacobson and MosesI (119841) proposed that the current diffu- 
sivity should be modified to include the diffusion of electrons 
across the large scale magnetic field due to the small scale 
field line stochasticity. The consequent modifications to the 
ohmic resistivity have only a marginal effect on the Sweet- 
Parke r reconnection speeds. iMatthaeus and LamkinI (119851 
1986h studied 2D magnetic reconnection in the presence of 
external turbulence both theoretically and numerically. They 
pointed out various turbulence mechanisms that would en- 
hance reconnection rates, including multiple X-points as re- 
connection sites. However, this work did not include the ef- 
fect of magnetic field wandering, which is at the core of the 
LV99. They did not provide analytical predictions of the re- 
connection speed eitheilj. 

We begin by offering a brief summary of the differ- 
ences betwe en the Sweet - Parker model of the laminar re- 
connection (IParken 1 19571: ISweetl Il958h and the Lazarian- 
Vishniac model which accoun ts for the effects of turbulence 
dLazarian and Vishniag Il999h . The latter can be seen as a 
generalization of the Sweet-Parker model (see Fig. [T) in the 
sense that the two regions of differing magnetic directions 
are pressed up against one another over a broad contact re- 
gion. This is a generic configuration, which should arise 
naturally whenever a magnetic field has a non-trivial con- 
figuration, whose energy could be lowered through recon- 
nection. The outflow of plasma and reconnected flux will 
fluctuate as the turbulence evolves and the field line connec- 
tions change, but the long term average will reflect the tur- 
bulent diffusion of the field lines. Consequently, the essen- 
tial difference between the Sweet-Parker model and the LV99 
model is that the former the outflow is limited by microphys- 
ical Ohmic diffusivity, while in the LV99 model the large- 
scale magnetic field wandering determines the thickness of 
outflow. For extremely weak turbulence, when the range of 
magnetic field wandering becomes smaller than the width of 
the Sweet-Parker layer LS~^^^, the two models are indistin- 
guishable. By weak turbulence, following LV99, we under- 



Sweet-Parker model 

L 



'A t the same time, after LV99 was published. iKim and Diamond 
1 120011) produced a study arguing that turbulence will not change re- 
connection reconnection rat es in the Sweet - Parke r geo metry. This 
study has been criticized by iLazarian et al.l (120041) and lEvink et al.l 
( 120111) . The present paper provides numerical evidence that the re- 
connection rates do increase in the presence of turbulence. 





Turbulent model 




Fig. 1. Upper plot: The Sweet-Parker reconnection model. The 
outflow is confined to a thin layer of 5, which is set by Ohmic diffu- 
sivity. The length of the cuiTent sheet is a macroscopic scale L^S. 
Magnetic field lines are assumed to be laminar. Middle plot: Re- 
connection in the presence of stochastic magnetic field lines. The 
stochasticity introduced by turbulence is weak and the mean field 
is clear direction. The outflow width is set by the diffusion of the 
magnetic field lines, which is a macroscopic process, independent 
of resistivity. Low plot: An individual small scale reconnection re- 
gion. The reconnection over small patches of magnetic field de- 
termines the local reconnection rate. The global reconnection rate 
is substantially larger as many independent patches reconnect si- 
multaneously. Conservatively, the LV99 model assumes that the 
small scale events ha ppe n at a slow Sweet-P arker rate. Following 
ILazarian et id] ( |2004|) and lKowal et"ai] J2009l) . 



Stand a regime where the correlation length is much greater 
than the distance by which individual field lines deviate from 
a straight line. 

LV99 considered a large scale, well-ordered magnetic 
field, of the kind that is normally used as a starting point 
for discussions of reconnection. In the presence of turbu- 
lence the field has some small scale 'wandering'. LV99 sug- 
gested that the presence of a random magnetic field compo- 
nent leads to fast reconnection. There are three phenomena 
mainly responsible for this: 

- only a small fraction of any magnetic field line is subject 
to direct Ohmic annihilation, therefore the fraction of 
magnetic energy that goes directly into heating the fluid 
drops down to zero as the fluid resistivity vanishes, 

- the presence of turbulence enables many magnetic field 
lines to enter the reconnection zone simultaneously. 
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- turbulence broadens the width of the ejection thickness 
allowing for more efficient removal of the reconnected 
flux. 



With the iGoldreich and Sridhaij d 19951 henceforth GS95) 
model of turbulence, LV99 obtained: 



Vrec = Va 



1/2 



VaJ 



(1) 



where I and vi are the energy injection scale and velocity. 
This expression assumes that energy is injected isotropically 
at the scale I smaller than the length of current sheet L, which 
for subAlfvenic turbulence leads to the generation of weakly 
interacting waves at that scale. The waves transfer energy 
to modes with larger values of k± until strong turbulence 
sets in. It is important to note that the strongly turbulent ed- 
dies have a characteristic velocity of vturb ~ Va{vi/Va)'^. In 
other words, the reconnection speed is the large eddy strong 
turbulent velocity times factors which depend on whether the 
current sheet is smaller or larger than the large eddies (whose 
length is approximately the injection scale). In this sense the 
reconnection speed should be fairly insensitive to the exact 
mechanism for turbulent power injection. The main purpose 
of this paper is test whether or not this is true for a simple 
modifi cation of the driving mechanism used in iKowal et al.l 
(l2009l) . 

It is important to note three features of Equation ([T). First, 
and most important, it is independent of resistivity. This is, 
by definition, fast reconnection. Second, we usually expect 
reconnection to be close to the turbulent eddy speed, the geo- 
metric ratios that enter the expression, i.e. the injection scale 
I divided by the length of the reconnection layer L, are typ- 
ically of order unity. Reconnection will occur on dynami- 
cal time scales. Finally, we note that in particular situations 
when turbulence is extremely weak the reconnection speed 
can be much slower than the Alfven speed. 

More recently. Equation ([T) was derived using the ideas 
based on the well- known concept of Richardson diffusion 
(jEyink et al., l201lb . From the theoretical perspective this 
new derivation avoids rather complex considerations of the 
cascade of reconne ction events that we re presented in LV99 
to justify the model. lEvink et al.l (1201 lb also shows that LV99 
model is closely connected to the recently developed idea of 
"spontaneous stochasticity" of magnetic fields in turbulent 
fluids. 

In general, the situation in reconnection community now is 
very different from that a decade ago. Currently, possibilities 
of fast reconnection in MHD regime due to instabilities of 
the reconnection layers are widel y discussed (ILoureiro et all 
2009t iBhattachariee e t al., 2009). These ideas can be traced 
back to the work of iShibata and Tanuma (2001). The in- 
stabilities, like tearing instability, open up the reconnection 
layer enabling a wide outflow. We expect such an outflow 
to become turbulent for most of astrophysical conditions. In 



this case, the process can be important for initiating recon- 
nection in the particular situation when the level of preex- 
isting turbulence is initially low to initiate sufficiently fast 
reconnection. We feel that exploring the ways of initiation of 
turbulent reconnection is very synergistic to the LV99 ideas, 
but in the current paper we focus on the case of preexisting 
turbulence of sufficient level. This is the primary domain for 
which LV99 provides predictions. 

Given the limited dynamical range of numerical simula- 
tions, we can only inject power on scales less than than L. 
The most convenient numerical parameter is not vi, but the 
energy injection power P. The power in the turbulent cas- 
cade is P '^ ''-'turbi^A/l) or vf/{lVA)- The amount of energy 
injected during one Alfven time unit Ia = L/Va, which is 
constant in our models, is i^P ~ {L /VA)vf / (IVa) ■ There- 
fore vf - {l/L)^^'^{PtAy^'^VA. Substituting vf in Equa- 
tion ([T]| results in 



\l/2, 



Vrec--{j^]{tAPr''oclP'^^ 



(2) 



which is the prediction we will test here. In what follows 
we refer to the injection power and scale using Pinj and Imj, 
respectively. 

3 Numerical Setup 

3.1 Governing Equations 

We use a high-order shock-capturing Godunov-type scheme 
based on the monotonicity preserving (MP) spatial re con- 
stmction (see Suresh and Huvnhlll997l:lHe et al .'. 2011, e.g.) 
and Strong Sta bility Preserving Run ge-Kutta (SSPRK) time 
integration ("see iGottlieb et al.[l2009l and references therein) 
to solve isothermal non-ideal MHD equations. 



dp 



-V-(pv)-0, 



dpv 



-V- 



pw- 



n2 \ 1 

a^p+ — / BB 

^ St: J 47r 

dA , 



=f, 



g, 



(3) 
(4) 
(5) 



where p and v are plasma density and velocity, respectively, 
A is the vector potential, E — ~v x B + j]J is the electric 
field, B = V X A is the magnetic field, 7 = V x B is the 
current density, a is the isothermal speed of sound, 77 is 
the resistivity coefficient, and / and g represent the turbu- 
lence driving terms either in velocity or vector potential. We 
used a multi-state Harten-Lax-van Leer (HELD, Mignong, 
2007h approximate Riemann solver for solving the isother- 
mal MHD equations. The HELD Riemann solver takes 
into account magnetic fields, and can follow Alfven waves 
with minimal numerical dissipation. This is particularly 
important here, because our simulations are in the quasi- 
incompressible regime, where most of energy is transported 
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by Alfven waves. The V JB = is maintained by solving the 
induction equation (Eq. |5]l using the field interpolated con- 
strai ned transport (CT) sch eme base d on a sta ggered mesh 
(e.g. iLondrillo and Del Zanna. .2000i: IfHflii I2OOQ) . 



3.2 Model Description and Initial Conditions 

Our reconnection simulation setup is illustrated in Figure |2] 
which is a 2D cut through the problem setup, indicating the 
location of the diffusion region. The top and bottom of the 
computational domain contain equal and opposite field com- 
ponents in the x direction, as well as a shared component 
Bz (see the left panel of Fig. |2]i. Magnetic field lines enter 
through the top and bottom and are bent by the inflow Vin 
as they move into the diffusion region. The diffusion region 
has a length A in the x direction and a thickness 6 in the y 
direction (see the left panel of Fig.|2|l. The box is periodic in 
the z direction and the diffusion region extends through the 
entire domain. The projection of the magnetic topology on 
the XZ plane shows that the lines in the upper region (solid 
lines in the right panel of Fig. |2]l and in the lower region 
(dashed lines) create an angle a determined by the strength 
of the shared component Bq^ . Once the incoming magnetic 
lines enter the diffusion region, they are reconnected and the 
product of this process is ejected along X direction with a 
speed Vout (the left panel of Fig.|2]i. 

We begin with a Harris current sheet of the form 
Bx{x,y,z) = i3o2;tanh(y/0) initialized using the magnetic 
vector potential Az{x,y,z) = ln|cosh(j//6')|,and a uniform 
guide field Bz{x,y,z) = Bq^ = const. The initial setup is 
completed by setting the density profile from the condi- 
tion of the uniform total (thermal plus magnetic) pressure 
Ptot{t = 0,x,y,z) = const and setting the initial velocity to 
zero. 

Magnetic reconnection is initiated by a small perturbation 
of the vector potential (5j4^(x,y,z) — SBoxCOs{2ttx) to the 
initial configuration of Az{t = 0,x,y,z) whose strength is 
given by the coefficient 5Bqx- 

In all our simulations the strength of the magnetic field is 
expressed in terms of the Alfven velocity defined by the an- 
tiparallel component of the unperturbed magnetic field. Sim- 
ilarly the density is expressed in terms of the unperturbed 
density po — l and velocities are expressed as fractions of the 
fiducial Alfven speed. The length of the box in the x direc- 
tion defines the unit of distance and time is measured in units 
of Lx/Va- In the new set of models we set the initial strength 
of the antiparallel magnetic field component Bq^ = 1.0 and 
the guide field B^z =0.1. We performed modeling for two 
resistivity coefficients riu — 5- lO^"* and rju = 10^"^ which are 
expressed in the dimensionless units. The initial perturbation 
is set to SBqx = 0.024. In order to avoid the complications 
of strong compressibility we have set the sound speed to 4.0. 



3.3 Boundary Conditions 

Our computational box has a grid of 256x512x256 or for 
higher resolution runs, 512x1024x512. In dimensionless 
units its size is L^ — Lz^l and Ly = 2. We double the size 
in the y direction to keep the driven turbulence away from 
the inflow boundary. There is no physical reason to do this, 
but driving turbulence near the inflow boundary produces nu- 
merical instabilities. 

As mentioned earlier, we use three different types of 
boundary conditions, depending on the direction of the 
boundary: outflow boundary conditions along the x direc- 
tion, inflow boundary conditions along the y direction and 
(sometimes) periodic boundary conditions along the z direc- 
tion. 

The open boundary conditions are the same as tho s e used 
in our previous modeling. We refer to iKowal et al.l (l2009h 
for their detailed description. Briefly, we use simple "zero- 
gradient" boundary conditions, setting the normal derivatives 
of the fluid variables (density and momentum) to zero. In the 
hydrodynamic limit this allows waves to leave the box with- 
out significant boundary reflections. In turbulence simula- 
tions this can lead to a slight drift in the fluid density. There is 
no requirement that the boundary density is constant, and in- 
flows and outflows can cause a small net gain or loss from the 
system. Fortunately, changes in the tota l mass are small an d 
only fluctuate around the initial value (IKowal et al.[ 120091) . 
They do not influence our results significantly. 

In order to incorporate the magnetic field into the open 
boundary conditions, we set the transverse components of 
the vector potential A using first order extrapolation. The 
normal derivative of the normal component is set to zero. 
In this way the normal derivatives of the transverse compo- 
nents of the magnetic field are zero, while the normal compo- 
nent of magnetic field is calculated from the zero-divergence 
condition V • B = 0. This approach avoids the generation of 
nonzero magnetic divergence at the boundary. However, it 
has the drawback that it creates a small jump in the momen- 
tum flux across the boundary resulting from the presence of 
non-zero terms {—Bx,By,Bz)dxBx at the X outflow bound- 
ary and {Bx,—By,Bz)dyBy at the Y inflow boundary. We 
have evaluated the velocity increment these terms produce at 
each time step. In models with the strongest turbulence these 
terms were of order of 10^^ and 10^^ at the X and Y bound- 
aries, respectively. In the presence of strong outflows and 
inflows, generally of order unity, they are clearly negligible. 

Simulations with explicit resistivity run into problems at 
the boundaries. In order to avoid a non-continuous resistive 
term and difficulties with the treatment of the current density 
/ we have introduced a zone of decaying resistivity near the 
boundary. In a thin layer near the boundary, the value of 
resistivity rju decays down to a very small value chosen to 
be close to the numerical resistivity rjn of our code. In our 
models we adopt the value of 77„ = 3 • lO^"'. None of this has 
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Vu^AAAAA/>AnA^U^A"AA/^AArt/WVWVVVl/^A'\AA^^ 



periodic 



Fig. 2. A schematic of our magnetic field configuration projected on the XY (left) and XZ (right) planes. Left: XY projection of the magnetic 
field lines. The gray area describes the diffusion region where the incoming field lines reconnect. The longitudinal and transverse scales of 
the diffusion region aregiven by A and 5, respectively. We use outflow and inflow boundary conditions in the x and y directions, respectively. 
Right: XZ projection of the magnetic field lines as seen from the top. Solid and dashed lines show the incoming field lines from the upper and 
lower parts of the domain, respectively. We see that the oppositely directed field lines are not antiparallel but are set an angle ct deter mined by 
the st rength of the shared component Bz. The z boundary conditions can be open or periodic, depending on the model. (from iKowal et all 
I2OO9I) 



an effect on the reconnect ion speeds. The v alidation of this 
method was presented in .Kowal et al.l (120091) . 



3.4 New Method of Turbulence Driving 

In our previ ous work we dr ove turbulence using a method 
described by lAlveliusI (1 1 999h . in which the driving term was 
implemented in the spectral space with discrete Fourier com- 
ponents concentrated around a wave vector kinj correspond- 
ing to the injection scale Unj ~ 1/kinj. We perturbed a num- 
ber Nf of discrete Fourier components of velocity in a shell 
extending from kinj ~ ^hnj to kinj + ^kinj with a Gaus- 
sian profile of the half width kc and the peak amplitude Vf 
at the injection scale. The amplitude of driving is solely de- 
termined by its power Pinj, the number of driven Fourier 
components and the time step of driving Atf. The random- 
ness in time makes the force neutral in the sense that it does 
not directly correlate with any of the time scales of the tur- 
bulent flow, and it also determines the power input solely by 
the force-force correlation. 

On the right hand side of Equation (|4|l, the forcing is repre- 
sented by a function/= pii, where p is local density and ii is 
random acceleration calculated using the method described 
above. In a similar way we can drive turbulence in the vector 
potential or magnetic field, which is represented by term g 
on the right hand side of the induction equation (Eq.|5]l. 

In the new method of turbulence driving we add individual 
eddies with random locations of their centers and random 
orientations, either to velocity or magnetic field, at random 
moments in time. This guarantees the randomness of new 
forcing. 

Each eddy is calculated from a kernel function described 
by a directional vector a (with amplitude \a\) multiplied by a 



Gaussian function 



^{r) =aexp I 



252 



(6) 



where re is the location of the eddy center, and S is the eddy 
width. An actual eddy is generated from such a kernel func- 
tion by taking its curl, i.e. Sf—pCVx ^) dt or Sg — V x'^dt 
in the case of injection in velocity or magnetic field, respec- 
tively. For example, if we assume that we inject one eddy 
in the magnetic field at re = (0,0,0), and that the perturbing 
vector potential fluctuation has only the nonzero component, 
i.e. $ = (0,0,$^), the contribution to magnetic field is ex- 
pressed by 



Sgy 



{x,y,z): 



<52 



exp 



252 



X 





dt . 



(7) 



This function describes an eddy injected in the XY plane with 
the maximum amplitude 6g„iax — \a\5^^e^2 at the distance 
Tmax = 5 and injection scale Unj = 5. We know the energy 
injected by one eddy, which is AEeddy = tt'^/^ \a\'^6/2, there- 
fore we can determine its amplitude \a\ from the injection 
power Pinj and the injection rate Ninj, which is the number 
of eddies injected in a time unit. 



Pinj 



■N^njAE, 



2P 



eddy 



171J 



7Ty^N,n,S 



(8) 



These estimates are done for the 3D case. In the 2D case the 
eddy energy is AEeddy = 7r|fl|^/2, and therefore the eddy 
amplitude can be determined from \a \ = \j2Pinj / {T^Ninj ) . 

In the new method there is no direct treatment of the 
velocity-force correlation, therefore there is no guarantee that 
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this correlation is zero and the injected power is completely 
determined by the force-force correlation. A reasonable so- 
lution to this problem would be to control the amount of in- 
jected energy and modify the amplitude of injected eddies or 
the injection rate at each time step in order to compensate 
the differences in the energy injected in the domain. The per- 
formed tests show, however, that although the velocity-force 
correlation is not zero, it is in fact fluctuating in time with a 
small amplitude, and giving in result a zero net contribution. 

The new method drives turbulence directly in the real 
space, in contrast to the previous one, therefore it can be 
applied locally. We drive turbulence in a subvolume of the 
domain. The size of the subvolume is determined by two 
scales, the radius r^ on the XZ plane around the center of 
the domain and the height hd describing the thickness of the 
driving region from the midplane. In this way we avoid driv- 
ing turbulence at the boundary and reduce the influence of 
driving on the inflow or outflow. 

All models are evolved without turbulence for several dy- 
namical times in order to allow the system to achieve station- 
ary laminar reconnection. Then, at a given time <{, we start 
driving turbulence by increasing its amplitude to the desired 
level, until te- In this way we let the system to adjust to a 
new state. From time te the turbulence is driven with the full 
power Pinj. 

3.5 Reconnection Rate Measure 

In the next sections we measure the reconnection rate using 
the new method of reconnection rate measure introduced in 
Kowal et alj ( 120091) and described by a formula 



K-c 



1 



'^\-'^X,00 1-^2 



<isign{B^)E-dl-dt f\B^\dA 



(9) 



where B^ is the strength of reconnecting magnetic compo- 
nent, E is the electric field, dA is area element of an XZ plane 
across which we perform integration, dl is the line element 
separating two regions of the YZ plane defined by the sign of 
Bx, l^oj.ool is the asymptotic absolute value of B^, and Lz 
is the width of the box. 

This method of the reconnection rate measure was derived 
from the magnetic flux conservation '5 and takes into ac- 
count all processes contributing to the change of magnetic 
flux. The electric field v xB — ry can be further divided into 
an advection term v x B^x, a shear term v x (Byij + Bzz), 
and a resistive term ~j]j. With this in mind the line integral 
can be rewritten as 



(10) 



(psign{Bx)E-dl— (p |-Ba;|(v^ xx)-dl 

+ (psign{Bx)vx{S:xBj_)-dl — (prj-dl. 



This new reconnection measure contains the time deriva- 
tive of the absolute value of B^, and a number of bound- 
ary terms, such as advection of B^ across the boundary and 
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Fig. 3. Evolution of total mass A/, and kinetic and magnetic ener- 
gies, Ekin and Emag, respectively. Two dotted vertical lines bound 
the period of gradually increasing turbulence. The resistivity in this 
model is set to r; = 10^'^ and the shared component of magnetic 
field Bqz ~ 0.1. In this model we inject turbulence in the magnetic 
field. 



the boundary integral of the resistive term rg. The additional 
terms include all processes contributing the time change of 
\Bx\- In particular, they can have non-zero values. 

3.6 Table of Simulated Models 

In Table [1] we list parameters of a ll the model s prese nted 
in this paper including models from lKowal et al.l (120091) and 
models with new driving. As in the previous paper we di- 
vided them into several groups. In each group we calculated 
models in order to study the dependence of the reconnection 
rate on a characteristic parameter of turbulence or resistiv- 
ity. We have studied the dependence of reconnection on the 
power of turbulence (models "PD"), injection scale (mod- 
els "SD"), and viscosity (models "VD"). Models with new 
driving are marker with a right arrow ( — >), and models with 
new driving and higher resolution are marked with a symbol 
"HR". 

Only the varying parameters are listed in the table, the 
strength of guide field Bqz, the uniform resistivity rju, the 
uniform viscosity i/^, the power of turbulence Pinj and its 



injection scale kinj, and the method of turbulence driving. 

All models presented in this section were calculated with 
the grid size Ax « 0.004 corresponding to the resolution 
256x512x256, except the model marked with symbol "HR", 
which was simulated with the resolution 512x1024x512 
(Ax « 0.002). 



4 Results 

4.1 Time Evolution of Energies and Reconnection Rate 

In Figure [3] we present an example of the evolution of to- 
tal mass, and kinetic and magnetic energies in a model with 
Pin — 1.0, kf — 8, and 77^ = 10^'^. We inject turbulence in 
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Table 1. List of models. 



Name Bn 



riu [10-3] v^ [10-3] 



Afcinj Driving Type 



PD 


0.1 


1.0 


0.0 


0.1 


8 


old in V 




0.1 


1.0 


0.0 


0.2 


8 


old in V 




0.1 


1.0 


0.0 


0.5 


8 


old in V 




0.1 


1.0 


0.0 


1.0 


8 


old in V 




0.1 


1.0 


0.0 


2.0 


8 


old in V 




1.0 


1.0 


0.0 


0.1 


8 


old in V 




1.0 


1.0 


0.0 


0.2 


8 


old in V 




1.0 


1.0 


0.0 


0.5 


8 


old in V 




1.0 


1.0 


0.0 


1.0 


8 


old in V 




1.0 


1.0 


0.0 


2.0 


8 


old in V 


— > 


0.1 


1.0 


0.0 


0.2 


8 


new in B 


— > 


0.1 


1.0 


0.0 


0.5 


8 


new in B 


— > 


0.1 


1.0 


0.0 


1.0 


8 


new in B 


HR 


0.1 


0.5 


0.0 


1.0 


8 


new in V 


SD 


0.1 


1.0 


0.0 


1.0 


5 


old in V 




0.1 


1.0 


0.0 


1.0 


8 


old in V 




0.1 


1.0 


0.0 


1.0 


12 


old in V 




0.1 


1.0 


0.0 


1.0 


16 


old in V 




0.1 


1.0 


0.0 


1.0 


25 


old in V 




1.0 


1.0 


0.0 


1.0 


5 


old in V 




1.0 


1.0 


0.0 


1.0 


8 


old in V 




1.0 


1.0 


0.0 


1.0 


12 


old in V 




1.0 


1.0 


0.0 


1.0 


16 


old in V 




1.0 


1.0 


0.0 


1.0 


25 


old in V 


— > 


0.1 


1.0 


0.0 


1.0 


8 


new in B 


— > 


0.1 


1.0 


0.0 


1.0 


24 


new in B 


— >• 


0.1 


1.0 


0.0 


1.0 


32 


new in B 


HR 


0.1 


0.5 


0.0 


1.0 


8 


new in V 


VD 


0.1 


1.0 


0.2 


1.0 


8 


old in V 




0.1 


1.0 


0.5 


1.0 


8 


old in V 




0.1 


1.0 


1.0 


1.0 


8 


old in V 




0.1 


1.0 


2.0 


1.0 


8 


old in V 




0.1 


1.0 


3.0 


1.0 


8 


old in V 




0.1 


1.0 


4.0 


1.0 


8 


old in V 




0.1 


1.0 


5.0 


1.0 


8 


old in V 



magnetic field using new forcing method, gradually increas- 
ing its strength from t — A to t — b. This period is marked 
by two dotted vertical lines in Figure |3] We see an increase 
of kinetic energy during this period due to the injection and 
saturation after i = 5. The kinetic energy preserves constant 
value during the turbulent stage very well. The magnetic en- 
ergy increases during this stage, slowly saturating. This in- 
crease is attributed to the injection of magnetic eddies. On 
the contrary, the total mass in the system decays slowly. We 
emphasize, that since we use open boundary conditions, not 
perfect conservation of mass and energies is possible in the 
presence of turbulence. 

In Figure |4] we show the evolution of reconnection rates 
Vrec for two modcls with the same set of initial conditions, 
but in the first model we drove turbulence by injecting mag- 
netic eddies using the new method described in this paper 



(black line), and in the second model we inj ect velocity 



fluctua tions using the old method described in iKowal et al 



(120091) (blue line). In this plot we recognize an increase 
of both rates during the introduction of turbulence. After 
the transition period between i = 4 and t = 5, when the sys- 
tem adjusts to a new state, both measures coincide and even 
though they fluctuate, they reach a stationary state character- 
ized by faster reconnection. Both types of turbulence bring 
the reconnection rate to similar level. Somewhat higher re- 
connection rate in the model with new driving could be at- 
tributed to the fact that this model was calculated using the 
5th order spatial reconstruction and the 3rd order integra- 
tion in time, in contrast to the old model where we used the 
second order methods. Lower order, especially in the spa- 
tial interpolation, introduces additional numerical diffusion 
decreasing the amplitudes of turbulent fluctuations at scales 
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Fig. 4. Evolution of tiie reconnection rate Vrec (black) for the same 
model as in Fig. [3] Blue line show the evolution of reconnection 
in a model with the same parameters in which the turbulence were 
driven using the old method. In this plot we present the measured 
rates of the Sweet-Parker reconnection VracSP and during the pres- 
ence of turbulence, Vrec,LV- Symbol SVrec,LV is the time variance. 
^Vrec,LV is the estimated uncertainty of the measure. 



comparable to the current sheet scale. 

In Figure|4]we also show the way of measuring the recon- 
nection rates, in the Sweet-Parker and LV99 stages, K-ecSP 
and Vrec,LV^ respectively. Because the reconnection rates 
fluctuate in the presence of turbulence we also measure their 
time variance dVrec,LV using the standard deviation. In ad- 
dition to the time variance of Vrec, we measure their errors 
by splitting the averaging region into two subregions and af- 
ter averaging the rates Virec and V2rec over each subregion 
(see Fig. |4]i, we take the absolute value of their difference 
AK-ec = Virec — V2rec- This difference corresponds to the 
error of Vrec, ie. it is different from zero if the rate is not 
constant in time. In all further analysis and presented plots 
we use values estimated in this way . These measures corre- 
spond exactly to those presented in .Kowal et alj ( 120091) . 



4.2 Topology of Magnetic Field 

In this section we compare field topologies in two example 
models run with the same set of parameters, but with dif- 
ferent types of driving. Both models have been simulated 
with the uniform resistivity rju = 10~^ and the resolution 
256x512x256. We injected turbulence with power Pj„j = 1.0 
at the injection scale kinj — 8. The only difference between 
models is the way we injected turbulence. In the old model 
we inject velocity fluctuations with random phases in Fourier 
space, and then transform them to real space and shape by a 
window in order to limit the injection to the specified region 
near the current sheet. In the new model, we inject magnetic 
loops with random locations and random orientations in the 
3D volume near the current sheet. The way of injecting tur- 
bulence is essentially different in both cases. 

In Figure|5]we show examples of XY-cuts (upper row) and 
XZ-cuts (lower row) through the box for the model with old 



driving. In the left and middle columns we show topologies 
of the velocity and magnetic field, respectively, with the in- 
tensities corresponding to the amplitude of components par- 
allel to the plotted plane. In the right column we show the 
absolute value of current density with overplotted magnetic 
vectors. Velocity has a very complex and mixed structure 
near the midplane due to constant injection of fluctuations in 
this region (see the left panel in Fig. |5]l. The majority of the 
velocity fluctuations is perpendicular to the mean magnetic 
field. This is because we are in the nearly incompressible 
regime of turbulence (large plasma /3) and most of the fluctu- 
ations propagate as Alfven waves along the mean magnetic 
field. Slow and fast waves, whose strengths are significantly 
reduced, are allowed to propagate in directions perpendicu- 
lar to the mean field as well. As a result, a big fraction of 
the turbulent kinetic energy leaves the box along magnetic 
lines. We observe, however, an efficient bending of magnetic 
lines at the midplane where the field is weaker (see the up- 
per middle plot in Fig. |5]l. This is not result of a driving, but 
result of reconnection. In general the interface between pos- 
itively and negatively directed magnetic lines is much more 
complex than in the case of Sweet-Parker reconnection. This 
complexity favors creation of enhanced current density re- 
gions, where the local reconnection works faster (see the 
right panel of Fig. |5]l. Since we observe multiple reconnec- 
tion events happening at the same time, the global reconnec- 
tion rate should be significantly enhanced. 

In Figure |6] we show similar examples of XY-cuts (up- 
per row) and XZ-cuts (lower row) through the box but for 
a model with the new way of driving turbulence. Here, a 
big number of individual eddies is injected in the magnetic 
field with random locations and random orientations in do- 
main. Comparing to plots in Figure |5] we see differences but 
also some clear similarities. Among the similarities, we note 
a highly turbulent region near the current sheet seen in all 
XY-cuts, with the current sheet itself strongly deformed and 
fragmented into many small scale current sheets (the right 
column of Figures |5] and |6ll. We see also some small in- 
crease of magnetic field strength near the current sheet (mid- 
dle top panels) resulting from working turbulence in the in- 
jection region. Among the differences we can list somewhat 
different distributions of the fragmented current sheets in the 
new model with clear enhancements in the locations where 
the magnetic eddies are injected at that moment. These en- 
hancements are clearly seen in the magnetic topology and 
current density plots (middle and right columns). In order 
to decrease those strong disturbances of the magnetic lines, 
we shall reproduce the same model with higher injection rate 
and reduced amplitudes of individual eddies. Another dif- 
ference is the strength of current density /. In the model 
with old driving we see more volume in which |7| reaches 
high magnitude and its structure is elongated with the local 
field. In the model with new driving, the current density with 
high strength seems to be less correlated with the local field, 
probably due to the presence of newly injected eddies. In the 
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Fig. 5. Topology and strength of the velocity field (left panel) and magnetic field (middle panel) in the presence of fully developed turbulence 
for an example model with old driving at time t = 12. In the right panel we show distribution of the absolute value of current density \J\ 
overlapped with the magnetic vectors. The images show the XY-cut (upper row) and XZ-cut (lower row) of the domain at the midplane of 
the computational box. Turbulence is injected with power Pi„j = 1 at scale ki-aj — 8 into velocity. Magnetic field reversals observed are due 
to magnetic reconnection rather than driving of turbulence, which is subAlfvenic. 



intermediate strengths, the structure of |/| seems to be better 
correlated with the local field. 

We see from this comparison that models with different 
driving of turbulence demonstrate different topologies of the 
fields. In the next sections we show, that the averaged recon- 
nection rates do not change significantly, confirming that the 
way we inject turbulence is of less importance and only its 
strength and injection scale have influence on Vrec- 

4.3 Dependence on Turbulence Strength 

Models with the new method of turbulence driving are listed 
in Table[T] We run a few models with varying turbulent pow- 
ers in order to verify if the new driving modifies our previous 
results. In these models we kept the same parameters as in 
the previous ones, which allowed to confirm the dependence 



of the reconnection rate Vrec on the power of injected turbu- 
lence P, 



tnj- 



Figure |7] shows the values of reconnection speed Vrec in 
models with turbulent power Pinj varying in the range of val- 
ues by more than one order of magnitude, from 0.1 to 2.0, for 
all previously shown models (black symbols) ini Kowalet al 



(2009) and for new models (blue and red symbols) in which 
we drove turbulence using the new method. Because the evo- 
lution of Vrec in new models reaches stationarity after time 
t = 6, we averaged Vrec from i = 6 to t = 10 in these models. 
Figure |4] shows that the reconnection rates oscillate around 
their mean values. In Figure |7] we plot how the averaged 
reconnection speed depends on the strength of turbulence. 
Filled symbols represent the averaged reconnection rate in 
the presence of turbulence. The dotted line corresponds to 
the reconnection rate during the Sweet-Parker process, i.e. 
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Fig. 6. Topology and strength of the velocity field (left panel) and magnetic field (middle panel) in the presence of fully developed turbulence 
for an example model with new driving at time t = 10. In the right panel we show distribution of the absolute value of current density \J\ 
overlapped with the magnetic vectors. The images show the XY-cut (upper row) and XZ-cut (lower row) of the domain at the midplane of 
the computational box. Turbulence is injected with power Pinj = 1 at scale ki„j = 8 directly in magnetic field. 



without turbulence. The error bars show the time variance 
of Vrec- The size of symbols indicates the uncertainty in our 
estimate of the reconnection speed AVrec,LV normalized to 
the uncertainty in the reconnection speed during the Sweet- 
Parker evolution AVrec.sP- It is calculated from a formula 
size = 2.0 — hiAVrec.Lv/^nAVrec,SP- ii AVrec.LV is of the 
order of AVrec,sp their symbols have the same sizes. 



The reconnection rates for models with new driving, which 
is described in ^3.41 confirm the theoretical dependence of 

1/2 

Vrec on the injected power, which scales as ~^ P^ ' . There 
is no significant difference between models in which turbu- 
lence was driven in velocity and in magnetic field. This is 
in agreement with the LV99 prediction, that the reconnection 
rate does not depend on the type of turbulence. 



4.4 Dependence on Injection Scale 

The reconnection rate Vrec in the presence of turbulence de- 
pends only on the strength of turbulence and its injection 
scale linj, according to Equation [T] for a fixed magnitude 
of the antiparallel magnetic field component. In the previ- 
ous subsection we presented studies on the turbulent power 
dependence. In this subsection we aim to study the injec- 
tion scale dependence. For this purpose we performed sev- 
eral models with the new way of driving turbulence, as well, 
in order to verify if they confirm the dependence of the re- 
connection speed Vrec on the scale linj at which we inject 
turbulence. The new models are listed in Table [T] We inject 
turbulence at several scales, from kinj = 8 to kinj — 32. At 
the upper end of this range the turbulence barely broadens 
the Sweet-Parker current sheet. At the lower end the turbu- 
lent eddies are barely contained within the volume in which 
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Fig. 7. The dependence of the reconnection speed Vrec on Pi„j 
updated by symbols from new models. Blue symbols show models 
with new driving in which the eddies where injected in magnetic 
field instead of velocity, as in the previous models (black symbols). 
The dotted line corresponds to the Sweet-Parker reconnection rate 
for models with ?)„ = 10"''. A unique red symbol shows the recon- 
nection rates from model with new driving in velocity performed 
with higher resolution (512x1024x512) and resistivity coefficient 
reduced to r^u = 5 ■ 10^*. Error bars represent the time variance of 
Vrec- The size of symbols corresponds to the error of Vrec (the way 
we calculate errors is described in ^4.1). 



we excite turbulent motions. 

In Figure |8]we present the reconnection speed dependence 
on the injection scale. We plot the averaged Vrec for old mod- 
els (black symbols) completed by the values from new mod- 
els with alternative driving (blue and red symbols). From the 
plot we clearly see a strong dependence of the reconnection 
rate on the injection scale. The new models very precisely 
follow the same dependence, confirming again that the type 
of turbulent driving has no influence on the process of re- 
connection, and only the power and injection scale of this 
driving have strong importance. Similarly, as in the power 
dependence plot, the new models have slightly higher recon- 
nection speeds comparing to the old ones. This is due to 
reduced numerical dissipation of velocity, since in the new 
models we used higher order methods. Dissipation removes 
energy at small scales. If it is smaller, due to higher order nu- 
merical scheme, the turbulent fluctuations reach higher am- 
plitudes at the current sheet scale. This influences the rate of 
individual reconnection events improving slightly the global 
reconnection rate Vrec- 

Figure [8] shows a bit weaker scaling with the injection 
scale than that predicted by LV99 model, i.e. Vrec ^ hnj- 
We see several possible sources for the discrepancy. For in- 
stance, the existence of a turbulent inverse cascade can mod- 
ify the effective Unj- In addition, reconnection can also mod- 
ify the characteristics of turbulence, such as the power spec- 
trum and anisotropy. We aim to study these problems in the 



Fig. 8. The dependence of the reconnection speed Vrec on Unj 
with additional models in which turbulence was driven in a new 
way, as described in 93.41 Similarly to Fig. |7] blue symbols show 
models with perturbed magnetic field, and red symbols correspond 
to a high resolution model with reduced uniform resistivity in which 
turbulence was driven in velocity. The dotted line corresponds to the 
Sweet-Parker reconnection rate for models with rju ~ 10"'^. Error 
bars and the size of symbols have the same meaning as in Fig.|7] 



future work. 



4.5 Dependence on Viscosity 



In iKowal et alj ( 120091) we performed studies of the recon- 
nection rate on the resistivity, both the uniform and anoma- 
lous ones, and we obtained great agreement with the Sweet- 

1/2 

Parker scaling Vrec ~ Vy- for the case without turbulence, 
and no dependence on resistivity in the presence of turbu- 
lence, as was predicted in LV99. In this section we performed 
additional studies of the reconnection rate dependence on 
viscosity. The dissipation scale of turbulent cascade is re- 
lated to the magnitude of viscosity. If the dissipation works 
at scales larger than the cuiTent sheet thickness, the turbu- 
lence cascade stops before reaching the current sheet and the 
global reconnection rate should be reduced. The reconnec- 
tion will be enhanced still by the broadened ejection region, 
allowing for more efficient removal of the reconnected mag- 
netic flux. 

In Figure |9] we show reconnection rates for models with 
varying viscosity coefficient. Although there is not pre- 
diction for this dependence in the LV99 model, we could 
test it numerically. In the Figure |9] we see a weak depen- 
dence Vrec ~ v~^^'^- This dependence might be also use- 
ful in understanding the reconnection speed differences be- 
tween models with the same set of parameters but different 
resolutions, or solved with different orders of the numeri- 
cal scheme. At low resolutions or low order schemes, the 
numerical viscosity is expected to be larger, thus we should 
observe reduced reconnection speeds in those cases. This is 
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Fig. 9. The dependence of the reconnection speed Vrec on the uni- 
form viscosity coefficient v. As explained in the text, the reconnec- 
tion speed is reduced with increasing value of v. The dotted line 
corresponds to the Sweet-Parker reconnection rate. Error bars and 
the size of symbols have the same meaning as in Fig.[71 



confirmed in Figures |7] and H] where we compare old models 
done with the second order scheme, and new models done 
with higher order schemes and higher resolutions. In those 
plots all new models demonstrate slightly higher reconnec- 
tion rates. 



5 Discussion 

5.1 LV99 in CoUisional and Collisionless Plasma 

The LV99 model was introduced for both collisional and col- 
lisionless media and it claimed that the microphysics of colli- 
sionless reconnection events does not change the resulting re- 
c onnection rat e s. Thi s point was subjected to further scrutiny 
inC 



Evink et alj ( 120111) who provided a thorough investigation 



of the problem and concluded that for most of astrophysi- 
cal collisionless plasmas the LV99 model should be applica- 
ble, provided that plasma is turbulent. With turbulence being 
ubiquitous in astrophysical conditions, this hardly constraints 
the applicability of the LV99 model. 

The LV99 model of reconnection is applicable to the colli- 
sional medium, such as the ISM, which is both turbulent and 
magn etized, and wher e the Hall-MHD reconnection does not 
work fYamadai 120071) . For instance, for Hall-MHD recon- 
nection to be applicable, it is required that the Sweet-Parker 
current sheet Ssp width is smaller than the ion inertial length 
di. Thus, the "reconnection criterion for media to be col- 
lisionless" is {L/diY/"^ /{(jJcTe) < 1' which presents a much 
severe constraint on the possible rate of collisions. As a result 
magnetic reconnection happens to be mediated by the Hall- 
MHD only if the extend of the contact region L (see Fig.[T]| 



does not exceed 10^^ cm. Magnetic fields in the ISM should 
interact over much larger scales. 

The LV99 model works in astrophysical environments to 
which the Hall-MHD reconnection is applicable, as well, like 
Solar corona, interplanetary medium, if the level of turbu- 
lence is high enough. The reconnection on microscales can 
happen fast, i.e. in the Hall-MHD fashion. This may not 
change, however, the global reconnection rate. The LV99 
model shows that even with relatively slow Sweet-Parker 
reconnection at microscales the global reconnection is lim- 
ited not by Ohmic resistivity, but by the rate of magnetic 
field wondering. We believe that the Hall-MHD local recon- 
nection of magnetic fields is taking place in the interplane- 
tary medium, which is being tested by local in-situ measure- 
ments, while the global reconnection rates are determined by 
magnetic field wandering as prescribed in LV99. 

5.2 Limitations of 2D Reconnection 

In the absence of quantitative model to be tested, simula- 
tions aimed at studying reconnection speed have been done 
in 2D, both for collisional and collisionless regimes. This 
allowed to achieve higher resolutions (compared to those 
contemporary available in 3D) but substantially constrained 
magnetic field dynami cs. For instance, the c l osest study 
to ours was done by Matthaeus and Lamkinl (Il985h (see 
also Matthaeu s and LamkinLll986l) . The authors studied 2D 
magnetic reconnection in the presence of external turbu- 
lence. An enhancement of the reconnection rate was re- 
ported, but the numerical setup precluded the calculation of 
a long term average re connection rate. A mor e recen t study 
al ong the approach in iMatthaeus and Lamkinl (119851) is one 
in lWatson et alJ (120071) . where the effects of small scale tur- 
bulence on 2D reconnection were studied and no significant 
effects of turbulence on reconnecti on were reported for the 
setup chosen by the authors. Later. IServidio et al.l (120101) re- 
done the modeling of 2D tur bulent reconnection following 
Matthaeus and Lamkinl (1 19851) with much higher resolutions. 
They used an advanced technique to detect precisely all X- 
points in the domain and then performed statistical studies 
confirming the Sweet-Parker relation for the reconnection 
rate as a function of X-point geometry. The development of 
different techniques to study magnetic reconnection is very 
important. Even though their modeling was limited to one 
type of highly superAlfvenic decaying turbulence (the initial 
uniform magnetic field was zero), they reported reconnec- 
tion rates with normalized values 0.1 — 0.3 and confirmed 
the importance of turbulence for modifying the character 
of magnetic reconnection and specifies heating and trans- 
port as the effect of particular significance, as well as for- 
mation of Petschek-type "X-points" in 2D turbulence. Due 
to the lack of large scale magnetic field configuration, their 
model represents a specific case, far from generic situation 
observed in the astrophysical objects where the mean and 
turbulent components of magnetic fields have comparable 
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strengths. Therefore, these studies ca nnot predict the globa l 
reconnection rate, as well. Moreover, IServidio et al.l (tZOlOh 
interpreted successful numerical confirmation of the LV99 



Numerical studies on Fermi acceleration in turbulent re- 



connection have a long history (e.g. iMatthaeus et al.L 11984 : 



model as a result of strong turbulence, although iKowal et al. 



(l2009h addressed this problem carefully showing that the am- 
plitudes of velocity fluctuations, both injected and obtained 
from spectra of developed and stationary turbulence, are frac- 
tions of Alfven speed. 

The fact that our study is in 3D is essential, as the LV99 
model is intrinsically three dimensional. The general pic- 
ture is of tangled field lines with reconnection taking place 
via a series of "Y-points" or modified Sweet-Parker sheets 
distributed in some fractal way throughout the turbulence. 
A large scale Sweet-Parker sheet will be replaced by a more 
fractured surface, but the current sheets will occupy a vani sh- 
ingly small fraction of the total volume and the field reversal 
will remain relatively well localized. The model predicts that 
the reconnection speed would be approximately equal to the 
strong turbulent velocity with a modest dependence on the 
ratio of the eddy length to the current sheet length. There 
should be no dependence on resistivity. The major results 
contained in our figures showing the dependence of the re- 
connection speed on resistivity, input power and input scale 
agree with the quantitative predictions of the LV99 model. 
We are not aware of any competing models to compare our 
simulations with. 

The major differences from the present study stem from 
the fact that we test a 3D model of reconnection, as the 
LV99 depends on effects, like field wandering, that hap- 
pen only in 3D. In order to show how different 2D and 3D 
worlds are, we performed similar studies to those presented 
in IKowal et al. (2009), but lim i ting th e do main to two dimen- 
sions (see iKulpa-Dvbel etall 1201 Ol) . In iKulpa-Dvbel et al 
( 120 101) we demonstrated that 2D magnetic reconnection in 
the presence of turbulence depends on the Ohmic resistivity, 
therefore, it is not fast. Also, the dependencies on the tur- 
bulent power and injection scales were significantly weaker 
than in the LV99. This dependence of 2D reconnection 
rate on Ohmic resistivity in the presence of turbulence, al 
though weaker than the Sweet-Parker relation Vr 



V- 



-1/2 



has been independently confirmed by iLoureiro et al.l (l2009b 
studies, performed with very different approach. These dif- 
ferences call for deliberation with simple extension of con- 
clusions coming from the 2D modeling to natural for mag- 
netic field fully three dimensional world. 

5.3 Applications of the LV99 Model 

Reconnection is one of the most fundamental processes in- 
volving magnetic fields in conducting fluids or plasmas. 
Therefore, the identification of a robust process responsible 
for reconnection has many astrophysically important conse- 
quences. Below we list a few selected implications of the 
successful validation of the LV99 model. 



Goldstein et al.','l986': Ambr osiano et all 1 19881: brake et al. , 
2006; .Hoshino. .2012) . In the Sweet-Parker model, it has 
been shown that particles can accelerate due to the in- 
duced electric field in the reconnection zone (ILitvinenko 
2003.) . This one — shot acceleration process, however, is 
constrained by the narrow thickness of the acceleration zone 
which has to be larger than the particle Larmor radius and by 
the strength of the magnetic field. Therefore, the efficiency 
of this process is rather limited. Besides, it also does not pre- 
dict a power-law spectrum, as generally observed for cosmic 
rays. Observations have always been suggestive that mag- 
netic reconnection can happen at a high speed in some cir- 
cumstances, in spite of the theoretical difficulties in explain- 
ing it. For instance, the phenomenon of solar flares suggests 
that magnetic reconnection should be first slow in order to 
ensure the accumulation of magnetic flux and then suddenly 
become fast in order to explain the observed fast release of 
energy. The LV99 model can naturally explain this and other 
observational manifestations of magnetic reconnection. Con- 
sider a particle entrained on a reconnected magnetic field line 
(see Figd}. This particle may bounce back and forth be- 
tween magnetic mirrors formed by oppositely directed mag- 
netic fluxes moving towards each other with the velocity 
Vrec- Each bounce will increase the energy of a particle in a 
way consistent with the requirements of the first-order Fermi 
procesjj (de Gouveia Dal Pino and Lazarian, 2001. 2003 : 



de Gouveia dal Pino and Lazarianl 120051; iLazariani uOOa) . 
This is in contrast to the second-order Fermi acceleration that 
is frequently discussed in terms of ac celerating particles b y 
turbulence generated by reconnection (ILa Rosa et al.[ 120061) . 
The numerical studies of the particle acceleration supporting 
these ideas have been already started (Kowal et al.. 201 lab) . 
An interesting property of this acceleration mechanism is 
that it is also potentially testable observationally. since the 
resulting spectrum of acc elerated particles is different from 



that ar is ing from a shock. | de Gouveia Dal Pmo arid Lazarian 
(1200 ih : Ide Gouveia dal Pino and Lazarianl (l2005 i) used this 



mechanism of particle acceleration to explain the syn- 
chrotron power-law spectrum arising from the flares of the 
microquasar GRS 1915H-105. 

Furt her apphcations can be found in the solar physics. Fol- 
lowing |ZweibeLaiid_Vamadal (120091) we note that solar flares 
inspired much of the earlier research on reconnection (see 



Pneumanl 119811 : iBastian et al.L Il998h . As the plasma in 



volved is substantially rarefied, the restrictive conditions for 



"Another way of understanding the acceleration of energetic 
particles in the reconnection process above is to take into account 
that the length of magnetic field lines is decreasing during recon- 
nection. As a result, the physical volume of the energetic particles 
entrained on the field lines is shrinking. Thus, due to Liouville's 
theorem, their momentum should increase to preserve the constancy 
of the phase volume. 
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the coUisionle ss reconnection are satisfied in this particular 
environment. ICassak et alJ (l2005h stated that bistable Hall 
reconnection can be important in this case. Stochastic re- 
connection provides an alternative explanation. Indeed, an 
important prediction of the LV99 model is related to the re- 
connection instability that arises in the situation when the 
initial structure of the flux prior to reconnection is laminar. 
Reconnection at the Sweet-Parker rate is negligible. This al- 
lows magnetic flux to accumulate. However, when the degree 
of stochasticity exceeds a threshold value, the reconnection 
itself should excite more turb ulence, creating a positive feed - 
back resulting in a flare (see Lazar ian and Vishniacl 120091) . 
The instability is a generic property of laminar field recon- 
nection in both coUisionless and collisional environments. 
Referring to the Sun, one may speculate that the difference 
between gradual and eruptive flares arises from the original 
state of magnetic field prior to the flare, at least in some spe- 
cific situations. In the case when the magnetic field is suffi- 
ciently turbulent the accumulation of magnetic flux does not 
happen and the flare is gradual. Similarly, the observed spa- 
tial spread of energy release during solar flares may be due to 
the spread of the region of turbulent fields once reconnection 
is initiated at one place. Recent observations demonstrate 
that gradual flares occur rather in regions with large scale 
and weak magnetic fields f or which Alfven times are large 
dShibata and Magarai l201lb . In light of that, the difference 
in Alfven times may explain different time scales in gradual 
and impulsive flares. Further research is necessary for estab- 
lishing the role of turbulence in changing the time scale of 
flare evolution. 

The LV99 model can find its application in the re- 
moval of magri e tic flu x from the star formation re- 
gions. IShu et al.l (1200 6) showed that magnetic field is re- 
moved from the star forming core cluster faster than it 
is allowed by the standard am bipo lar diffusion sce nario 
dTassis and Mouschoviasl l2005allbl) . IShu et alJ (l2007h pro- 
posed a mechanis m using efficient r econn ection through 
"hyper-resistivity". ISantos-Lima et al.l (120101) performed nu- 
merical studies of such a concept, replacing the "hyper- 
resistivity" with efficient stochastic reconnection. They re- 
ported removal of strong anticorrelations of magnetic field 
through "reconnection diffusion", which can mimic the ef- 
fect of enhanced Ohmic resistivity. 

LV99 showed that fast reconnection of stochastic mag- 
netic field makes the models of strong MHD turbulence self- 
consistent, because the critical balance in the GS95 model 
requires the existence of eddy-type motions perpendicular to 
the magnetic field. In the absence of reconnection this would 
result in unresolved knots that should drain energy from the 
cascade. The estimates in LV99 showed that the rates of re- 
connection predicted by the model are sufficient to resolve 
magnetic knots within one period. 



5.4 Reasons for Slow Adaptation of the LV99 Model 

The LV99 model of magnetic reconnection in the pres- 
ence of weakly stocha sti c mag netic fields was proposed by 
Lazarian and Vishniad in 1 19991 However, due to a few ob- 
jective factors it met less enthusiasm in the community than, 
for example, the X-point coUisionless reconnection. We be- 
Ueve that there were three major factors responsible for this. 

1 . The coUisionless X-point reconnection was initiated and 
supported by numerical simulations, while LV99 was a 
theory. Its numerical testing became possible only re- 
cently. The reconnection subject had a history of failed 
theories and models, which without direct numerical 
support were taken with a grain of salt. 

2. The acceptance of the idea of astrophysical fluids gener- 
ically being in turbulent state had much less observa- 
tional support at that time compared to the present day. 
By now we have much more evidence which allows to 
claim that models not taking the pre-existent turbulence 
has little relevance to astrophysics. 

3. The analytical solutions of LV99 were based on the use 
of GS95 model of turbulence. The GS95 model of tur- 
bulence, in fact, was extended LV99 by introducing the 
concept of local reference frame for turbulent eddies 
and by extending the GS95 scalings to the subAlfvenic 
case. The GS95 theory was far from being generally 
accepted at the time of the LV99 publishing. 

The situation has changed substantially by now. 
First of all, GS9 5 was su ccessfully tested num e rically 



dCho and Vishniad I2OOOI: iMaron and Goldreichl I2OOI ; 



Cho et al.L 12002b and their ideas have 
to describing the Alf venic cascade in 
MHD turbulence (see ICho and LazarianL 



(see [^ 



been extended 
comp r essible 



2002, 2003 



Kowal and Lazarian, 2010u The so-called "Big Power Law 
in the Sky" indicating the presence of turbulence on scales 
from tens of parsecs to thousands of kilo meters has been 
extended (IChepurnov and Lazarian! I2OIOI) . and the obser- 
vations of gas and synchrotron emission provided extended 
number of direct turbulence measurements confirming 



their presence 



turbulence 
(see Burkh art et al. , 



20101 iGaensler et al. 



2011b . Finally, the situation has changed with the numerical 
testing of the LV9 9 model. The 3D MHD simulations in 



Kowal et alJ ( 120091) supported the predictions in the LV99 



paper and our present work goes further with testing this 
model, by including different types of energy injection. 



There are attempts to modify GS95 theory by supplementing it 
with additional effects, like dynamical alignment (^Boldvrevl uOOSi , 
u006), polarization (Beresnyak and Lazarian, 2006), non-locality 
iGogoberidzd ( l2007h . All these attempts, however, do not change 
the very nature of the GS95 model. Moreover, some rece nt studies 
iBeresnvak and LazarianI ( l2009l l2010h ; iBeresnvakI (l201lh indicate 
that the numerical motivation for introducing these attempts may be 
due to the insufficient inertial range of the simulations involved. 
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It is worth noting also, that there is some implicit obser- 
vational evidence in the favor of the LV99 model, like ob- 
servations of the thick reco nnection current ou tflow regions 
observed in the So lar flares dCiaravella and Ravmond, .20081 ) . 
Sych et alJ (l2009h . explaining quasi-periodic pulsations in 
observed flaring energy releases at an active region above 
the sunspot, proposed that the wave packets arising from the 
sunspots can trigger such pulsations. They established a phe- 
nomenological relation between oscillations in a sunspot and 
pulsations in flaring energy releases. This phenomenon can 
be naturally explained by the LV99 model. 



6 Conclusions 

In this article we performed additional testing of the LV99 
model of fast reconnection under different types of turbulent 
driving using 3D numerical simulations. We have introduced 
a new method of driving turbulence by direct injection of the 
velocity or magnetic eddies with random locations in the do- 
main. We analyzed the dependence of the reconnection speed 
on the turbulence injection power, on the injection scale, as 
well as on the viscosity. We found that: 

- We observe similar changes of the topology of the mag- 
netic field near the interface of oppositely directed mag- 
netic field lines in models with two different turbulence 
injection mechanisms. These changes include the frag- 
mentation of the current sheet, favoring multiple simul- 
taneous reconnection events, as well as a substantial in- 
crease in the thickness of the outflow of reconnected 
magnetic flux and matter. 

- The relation between the reconnection rate Vrec and tur- 
bulent power Pinj remains unchanged under two differ- 
ent mechanisms of energy injection and is confirmed by 

new models to be Vrec ~ Pinj ~ Vi^, in agreement with 
the LV99 prediction. Moreover, the injection in mag- 
netic field produces similar effects on the reconnection 
as injection in velocity, remaining the dependence unal- 
tered. 

- The reconnection rate grows with the size of the in- 
jected eddies, which can be directly related to the turbu- 
lence injection scale. The rate of growth, for the models 
with old and new driving mechanism, is approximated 
by Vrec ^ linj Scaling which agrees with the previously 
obtained scaling. Somewhat steeper LV99 prediction. 



Vre 



linj, could results from limitations in the dy- 



namic range available for study. 

Reconnection in the presence of weak turbulence is only 
weakly sensitive to viscosity i'. From performed sim- 
ulations we obtained a dependence Vrec '^ v^^^^ for 
one set of parameters: Pinj = 1.0, hnj =8, rju — 10"'^, 
Bo, =0.1. 
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